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ABSTRACT 

Objectives: Although well established, the association 
between socioeconomic position and health and health 
behaviour is not clearly understood, and it has been 
speculated that familial factors, for example, 
dispositional factors or exposures in the rearing 
environment, may be underlying the association. The 
objective was to compare prescription fillings within 
twin pairs who are partly or fully genetically identical 
and share childhood exposures. 
Design: Twin cohort study. 
Setting: Denmark. 

Participants: Data from the Danish Twin Registry 
were linked to registers in Statistics Denmark and the 
Danish Registry of Medicinal Product statistics. A total 
of 8582 monozygotic (MZ) and 15 788 dizygotic same 
sex (DZSS) twins were included. 
Outcome measures: Number of prescription fillings 
during follow-up (1995-2005) was analysed according 
to education and income. Results of unpaired and 
intrapair analyses were compared. 
Results: An inverse social gradient in filling of 
prescriptions for all-purpose and system-specific drugs 
was observed in the unpaired analyses. In the intrapair 
analyses, associations were attenuated some in DZSS 
and more in MZ twins. Filling of drugs targeting the 
nervous system was still strongly associated with 
income in the intrapair analyses. 
Conclusions: Familial factors seem to account for 
part of the observed social inequality in filling of 
prescription medicine. 
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INTRODUCTION 

Social inequalities in health and health behav- 
iour have been widely demonstrated in many 
societies, including the Scandinavian coun- 
tries. 1 Often social disparities in health are not 
confined to the most marginalised groups of 
society but are expressed as a gradient over the 
entire spectrum of social stratification. 2 The 
inverse relationship between socioeconomic 
position (SEP) and health has consistendy 
been observed across different social 



Strengths and limitations of this study 



■ This study is based on nation-wide register data 
covering more than a decade, and we present 
results from a largely unselected population of 
more than 8000 monozygotic and 15 000 dizyg- 
otic same sex twins. 

■ This is the first twin study on health inequalities 
to report such clear-cut findings with sufficient 
power for meaningful interpretation of intrapair 
analyses and for comparison of zygosity-specific 
results. 

■ Prescription fillings are not a perfect measure of 
health status, but we still argue that it is an 
inventive and fairly valid indicator in a Danish 
setting where prescription medicine has to be 
prescribed by a medical doctor and accounts for 
90% of the total drug sale. 

■ The twin study is clearly a powerful approach, 
but when interpreting findings, alternative expla- 
nations for the pattern of attenuation should be 
considered, including measurement error and 
unshared confounding. However, in this case we 
find it unlikely that these potential sources of 
bias could fully explain our findings. 



indicators and for numerous types of health 
outcomes. 3 4 Use of prescription medicine is 
one health-related outcome that has been 
found to be socially patterned, where Danish 
data have shown a greater use among people 
with low SEP compared with people with high 
SEP.'' b Denmark has a tax-financed decentra- 
lised healthcare system with a partly need- 
dependent reimbursement system for out- 
patient prescription drugs. All prescription 
drugs have to be prescribed by a medical 
doctor and, compared with a number of other 
countries, prescription medicine constitute 
the majority of the total drug sales in 
Denmark (approximately 90%). 7 Many studies 
have shown substantially reduced social differ- 
ences in drug use when health status have 
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been taken into account,' which may indicate that 
filling of prescription medicine to a large extent reflects 
health status. This is supported by the fact that comorbid- 
ity indices, based on prescription drug dispensings, have 
been established as strong predictors of mortality 10 and 
that drug use has proven to be a valid indicator of self- 
rated health. 11 12 However, other factors such as health- 
care seeking behaviour and access to healthcare 13 14 may 
also influence drug use, but the social gradient in these 
factors would, theoretically, be expected to be opposite 
that of health status, that is, greater drug use among 
people with high SEP compared with low SEP. Such a posi- 
tive association is generally not supported by the existing 
Danish literature. 5 6 1 

Although well established, the inverse association 
between SEP and health and health behaviour is not 
clearly understood, but many factors during the life 
course may play a role in the production of these social 
health differences. Suggested pathways include: material 
deprivation, behavioural factors and psychological 
factors. 16 It has also been suggested that it is not SEP in 
adulthood per se that influences health, but underlying 
familial factors already at play in early life, 17 18 for 
example, rearing environment or genetic make-up are 
associated both with the obtainment of adult SEP and 
health status later in life. 

In a recent commentary, Oilman and Loucks 19 discuss 
why childhood environment may be such an important 
underlying factor in the production of social health dif- 
ferences. They argue that because childhood encom- 
passes developmentally sensitive periods for the 
acquisition of cognitive, psychological and social skills 
that have a profound impact on an individual's capacity 
for long-term health and provide opportunities (and 
constraints) for adult socioeconomic conditions, it can 
potentially influence the association between any adult 
risk factor and later life health. 19 The potential role of 
genetic setup in explaining health inequalities is 
reviewed by Mackenbach, 20 who suggests that the most 
plausible set of hypotheses relates to genetic influences 
of personal attributes such as cognitive ability, personal- 
ity and bodily and mental fitness that may influence 
both SEP and adult health status. 

With the discordant twin pair design, it is possible to 
investigate the effect of adult SEP on prescription fillings 
eliminating confounding from factors shared by a pair of 
twins, exploiting that twins are pardy or fully genetically 
identical and have experienced a similar childhood envir- 
onment due to their common upbringing. This approach 
has been applied in a number of previous twin studies 
investigating different health outcomes, including mor- 

a i go go 

tality, hospitalisations, breast cancer and also softer 
endpoints such as depression scores and self-rated 
health. 24 25 However, the majority of these studies suffer 
from power limitations either due to small sample sizes 
and/or rare outcomes, making interpretations difficult. 

In this study we aimed to investigate if there was an 
effect of SEP in adulthood on the number of 
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prescription fillings above and beyond the effect of 
shared familial factors. Filling of prescription medicine 
is a quantitative health measure rendering sufficient 
power for the intrapair comparisons. Differences in 
genetic relatedness between monozygotic (MZ; genetic- 
ally identical) and dizygotic (DZ) twins (share on 
average 50% of segregating genes) may further give an 
indication of the possible types of confounding (ie, 
genetic or environmental factors). In the theoretical 
situation of no measurement error and no confounding 
from unshared factors, genetic confounding would be 
indicated if an observed social gradient in prescription 
medicine was pardy attenuated in DZ twins and fully 
attenuated in MZ twins, when shared familial factors had 
been taken into account. Likewise, an attenuation of 
similar size in DZ and MZ twins would be compatible 
with shared environmental confounding. Finally, to 
support a SEP effect on prescription fillings not due to 
confounding from shared familial factors, the associ- 
ation would have to persist in both DZ and MZ twins. 

To support the utility of prescription medicine as a 
health indicator we also analysed its association with 
mortality. 

METHODS 
Study population 

Data from the Danish Twin Registry were linked with 
information in the Danish Registry of Medicinal 
Products and administrative registers in Statistics 
Denmark. The study population consisted of 24 370 
twins, including 8582 MZ and 15 788 dizygotic same sex 
(DZSS) twins, who were born during 1921-1965 and 
alive and resident in Denmark at the beginning of 
follow-up (1 January 1995), after exclusions due to tri- 
plets (n=384) unknown zygosity (UZ; n=3731), and 
dizygotic opposite sex pairs (n=16 710), missing data on 
education or income (n=845) and non-intact twin pairs 
in which one twin was dead, emigrated or had missing 
information (n=2847). To evaluate the representative- 
ness of the results from the unpaired twin analyses, we 
also analysed a 5% population-based sample of the 
Danish population (n=137 300). 

Social indicators 

We used information on the highest attained education 
in 1995 at age of 30+ of the study population, cate- 
gorised according to the International Standard 
Classification of Education of 1997 26 into: primary/ 
lower secondary education (basic compulsory educa- 
tion), upper secondary education (secondary) and ter- 
tiary education (Bachelor's degree and above). Income 
was measured as an average of the equivalised gross 
household income (Danish currency) in 1994-1995, 
grouped into quartiles within each birth cohort. 
'Equivalised' income is adjusted to take into account 
that households can share resources. Inspired by the 
standard Organisation for Economic Co-operation and 
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Development equivalence scale, the household- 
denominator was constructed according to an equiva- 
lence scale that assigned a value of 1 to the first house- 
hold member, of 0.7 to the second household member 
and of 0.5 to each additional household member. 
Calculating quartiles by birth cohort is aiming at taking 
into account that income is not fully comparable at dif- 
ferent ages during the life course. An implication of this 
is that the income measure reflects a relative rather than 
absolute income level. 

Prescription medicine 

The Danish Registry of Medicinal Products was estab- 
lished in 1995 and contains information on the total 
sale of prescription medicine in Denmark registered by 
the pharmacies via a computerised accounting system. 
The drug information is recorded according to the 
Anatomical Therapeutic Chemical classification system 
(ATC) 28 where the active substances in drugs are divided 
into different groups, primarily according to the organ 
system on which they act. Other information includes 
date of dispensing and daily defined doses (DDD) , defined 
as the assumed average maintenance dose per day for a 
drug used for its main indication in adults. 28 Outcomes 
of interest were defined as the number of prescription 
fillings during follow-up for all-purpose medication 
(ACT-groups: A, B, C, D, G, H, J, L, M, N, P, Q, R, S, V, 
excluding Q which is for veterinary use, and G02B and 
G03A, which are for contraceptive use) and four groups 
of system-specific medication: ATC-C (heart/ circulatory 
system), ATC-N (nervous system), ATC-R (respiratory 
system) and ATC-A (alimentary tract/metabolism) . If no 
ATC code was assigned a prescription, it was excluded 
(1%). DDD was used as a supplementary outcome 
measure. 

Other covariates 

Sex and age were included in the analyses. Analyses 
were performed separately for MZ and DZSS twins. 
Zygosity had previously been determined by question- 
naire, 29 a method that has proven valid with an overall 
misclassification frequency of less than 5% validated 
against a classification based on genetic markers. 30 

Data analysis 

Analyses were carried out using Cox regression for recur- 
rent events ( proportional means regression) , thereby 
accommodating multiple events (prescription fillings) 
per individual. 31 Age was used as underlying time- 
variable and participants were followed from age at base- 
line to the age at the end of follow-up (31 December 
2005) or death or emigration, and the timing of each 
prescription filling during this period was recorded. 
Both unpaired and intrapair analyses were performed. 
In the unpaired analyses robust SEs were used to 
account for the interdependence of observations within 
pairs and within each individual. Results of the unpaired 
analyses were compared with the intrapair analyses, 



which were performed using a stratified Cox model: 
ENij(t)=^oj(t)e p ' ui . Here ENjj(t) denotes the mean 
number of prescription fillings at time t for an individ- 
ual i belonging to twin pair j, (J.qj(t) the pair-specific 
baseline mean, and e 13 the effect of the covariates 
adjusted for shared familial factors. Thus, in this analysis 
the mean number of prescription fillings was compared 
within twin pairs, thereby controlling for familial factors 
shared by a pair of twins. In the stratified model only 
SEP-discordant twin pairs contribute information to the 
fJ-parameter for SEP, although concordant pairs may 
contribute information to other (3-parameters. 

As counting process models do not allow simultaneous 
events, only one prescription of interest per day was 
included in the analysis. 

Analyses were performed separately for zygosity and 
sex, and potential interactions between the social indica- 
tors and zygosity, sex and age, respectively, were evalu- 
ated by including interaction terms in the models and 
using Wald test statistics. To analyse the association 
between prescriptions and mortality, we used a Cox 
regression model, continuously updating the number of 
prescription fillings over time. Descriptive intrapair cor- 
relations were calculated as polychoric correlation 
coefficients. 

Subanalyses 

Since DDD may reflect the actual drug use more pre- 
cisely than the number of prescription fillings, mean 
DDD were also analysed by means of a fixed-effects 
linear regression model taking censoring into account. 32 
This was treated as a subanalysis due to the highly right- 
skewed distribution of DDD, making linear regression a 
less optimal choice of model. 

Permission for linkage and use of data for this study 
was obtained from the Danish Data Protection Board 
(2000-54-0047). 

RESULTS 

According to table 1, showing the distribution of educa- 
tion and income in the population, men were better 
educated than women. Generally, the education of the 
twin population was comparable to that of the 
population-based sample. Only UZ twins had a markedly 
lower educational status and a higher proportion of 
missing information. The table also shows the mean 
income within each quartile. Twins have a somewhat 
higher mean income as compared to the population- 
based sample. 

The proportion of education-discordant twin pairs was 
30% for MZ twins and 42% for DZSS twins, while 61% 
MZ twins and 68% DZSS twins were discordant on 
income (table 2). 

In the twin population, 1.9% did not redeem any drug 
prescriptions during the total follow-up period and 25% 
redeemed less than one prescription/year, on average 
(results not shown) . Overall, women filled more 
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Table 1 Descriptives on educational status and income in a 5% population-based sample (n=144 749) and a population of Danish twins (n=48 887), including intact and 
non-intact twin pairs, stratified on sex and zygosity 





5% Population sample 




All twins 






DZSS 




MZ 




DZOS 




UZ 






N 


Per cent* 


Mean 


SD 


(n=48 887) 


Mean 


SD 


(n=18 426) 


(n=9636) 


(n=16 710) 


(n=3731) 


Sex 










N 


Per cent * 






N 


Per cent * 


N 


Per cent * 


N 


Per cent * 


N 


Per cent * 


Educational statusf 
































Men 


































Primary 


iLO SOU 


32 






8909 


35 






3553 


36 


1517 


31 


2882 


34 


957 


43 


Secondary 




43 






11 237 


44 






4309 


44 


2292 


47 


3834 


45 


802 


36 


Tertiary 


1 Q TOO 


19 






4515 


18 






1685 


17 


963 


20 


1530 


18 


337 


15 


Missing 


4UY iL 


6 






800 


3 






310 


3 


128 


3 


208 


2 


154 


7 


i oiai 


vo pon 


100 






25 461 


100 






9857 


100 


4900 


101 


8454 


99 


2250 


101 


Women 


































Primary 


Q.1 01 K 


43 






10 182 


44 






3857 


45 


1903 


40 


3633 


44 


789 


53 


Secondary 


oo onn 

<iO Buy 


33 






7805 


34 






2868 


34 


1730 


37 


2807 


34 


400 


27 


Tertiary 


h /i oon 


20 






4613 


20 






1696 


20 


1004 


21 


1705 


21 


208 


14 


Missing 


OD1 P 


4 






442 


2 






148 


2 


99 


2 


111 


1 


84 


6 


Total 


72 069 


100 






23 042 


100 






8569 


101 


4736 


100 


8256 


100 


1481 


100 


Income* (DKK) 


































Men 


































1st quartile 


17 348 


24 


88 537 


OD 0"7Q 


5931 


24 


QO oon 


07 CCO 


2251 


23 


1061 


22 


1929 


23 


690 


31 


2nd quartile 


17510 


24 


145 676 


0£ BQ7 


6523 


24 


1 AO Q£"7 
I4y OD/ 


OA 1 7Q 
tL<\ I / y 


2563 


26 


1208 


25 


2201 


26 


551 


24 


3rd quartile 


17 976 


25 


191 360 


34 863 


6769 


25 


1 QR 7SR 


31 365 


2612 


27 


1375 


28 


2246 


27 


536 


24 


4th quartile 


18 228 


25 


31 1 232 


141 784 


5960 


25 


314 180 


225 459 


2340 


24 


1218 


25 


2006 


24 


396 




Missing 


1618 


2 






278 


2 






91 


1 


38 


1 


72 


1 


77 




Total 


72 680 


100 






25 461 


100 






9875 


101 


4900 


101 


8454 


101 


2250 


100 


Women 


















2174 


25 


1158 


24 


2200 


27 


538 


36 


1st quartile 


18 958 


26 


88 201 


26 316 


6070 


26 


91 205 


25 698 


2316 


27 


1222 


26 


2186 


27 


353 


24 


2nd quartile 


18219 


25 


142 746 


27 689 


6077 


25 


146 427 


25 562 


2152 


25 


1224 


26 


2106 


26 


313 




3rd quartile 


17 753 


25 


18 942 


36 067 


5795 


25 


194 968 


33 291 


1882 


22 


1097 


23 


1730 


21 


232 


\ 


4th quartile 


16 126 


22 


308 876 


225 888 


4941 


22 


312 063 


173 113 


45 


1 


35 


1 


34 


0 


45 




Missing 


1013 


1 






159 


1 






8569 


100 


4736 


100 


8256 


101 


1481 


100 


Total 


72 069 


99 






23 042 


99 






2251 


23 


1061 


22 


1929 


23 


690 


31 



*May not sum to 100 because of rounding. 
-(-International Standard Classification of Education. 
^Average of equivalised gross household income in 
DKK, Danish crowns; DZOS, dizygotic opposite sex 



1994-1995 in quartiles. 

twins; DZSS, dizygotic same sex twins; MZ, monozygotic twins; UZ, unknown zygosity. 
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Table 2 Descriptives on intrapair discordances on educational status and income, stratified on zygosity in a population of 


Ranich t\A/inc (ri—OA ^70\ 












DZSS 




MZ 




(n=15 788) 




(n=8582) 






N 


Per cent* 


N 


Per cent* 


Educational status 










Primary vs secondary 


ooyu 


iLo 


1 4oU 


1 4 


Primary vs tertiary 


yuu 


c 
O 


ol D 




Secondary vs tertiary 
Any discordance 


OA Q/l 
t>004 


1 A 
AO 


^oyo 


10 
10 

on 
oU 


Income 










ist quamie vs tina qudnnt; 




l O 


I U / eL 


12 


1st quartile vs 3rd quartile 


1706 


11 


754 




1st quartile vs 4th quartile 


1158 


7 


490 




2nd quartile vs 3rd quartile 


2266 


14 


1194 


i 


2nd quartile vs 4th quartile 


1574 


10 


682 




3rd quartile vs 4th quartile 


2000 


13 


1080 




Any discordance 


10 752 


68 


5272 


61 


*May not sum to 100 because of rounding. 



prescriptions than men for all age groups and the 
number of prescription fillings increased with age in 
both sexes (table 3). 

ATC-N drug prescriptions were the most frequently 
redeemed in the youngest age group. The mean 
number of prescription fillings in the twin population 
was similar to that in the population-based sample, 
although slightly higher in the population-based sample. 

In order to get a better idea of the clinical outcomes 
underlying the observed prescription patterns, we have 
included a more detailed description of the most inci- 
dent types of drugs at the 3-digit and 4-digit ATG level 
(see online supplementary table SI). From the table it 
appears that the ATC-C category primarily comprises 
drugs indicated for hypertension and ischaemic heart 
disease, while the ATC-N category includes drugs indi- 
cated for a wide range of different diagnoses, including 
chronic pain disorders, insomnia, anxiety, psychoses, 
depression, etc. 28 In the ATC-R category the main types 
of drugs are drugs for obstructive airway disease such as 
asthma and chronic obstructive lung disorder 
(COPD). 28 Finally, drugs for acid-related disorders and 
diabetes are the main constituents of the ATC-A cat- 
egory. The results of an analysis of the association 
between SEP and prescription fillings for all-purpose 
and system-specific drugs are presented in table 4. 

For all-purpose drugs, MZ twins in the highest income 
quartile had a 30% lower mean number of prescription 
fillings than MZ twins in the lowest income quartile 
(MR=0.70 (95% CI 0.65 to 0.77)). Similarly, MZ twins 
with a tertiary education had 17% lower mean number 
of prescription fillings than twins with a primary educa- 
tion (MR=0.82 (95% CI 0.75 to 0.89)). The association 
was slightly less strong for ATC-C drugs while stronger 
associations could be observed for the other system- 
specific drugs, particularly for ATC-N. In the unpaired 



analyses, no large differences could be observed 
between MZ and DZSS twins and there was no strong 
evidence for an interaction between SEP and zygosity in 
the unpaired analysis (education p=0.81, income 
p=0.95) . In contrast, the effect of SEP was moderated by 
zygosity in the intrapair analyses (education p=0.02, 
income p=0.05), suggesting a more attenuated effect of 
SEP in MZ twins than in DZSS twins. In MZ twins, gener- 
ally no statistically significant effect of education could 
be observed, while income was still somewhat associated 
with all-purpose medication and quite strongly with 
ACT-N drugs. In DZSS twins, a minor or moderate 
attenuation could be observed, and all effect parameters 
were statistically significant, except for ATC-C drugs. In 
the population-based sample, the associations were 
similar to that in the overall twin population. 

When using DDD (for all-purpose drugs) as outcome 
measure (table 5), results similar to the results based on 
the number of prescription fillings were obtained, 
although the attenuation patterns in MZ twins were 
somewhat more clear in this analysis, showing no statis- 
tically significant associations with DDD of neither edu- 
cation nor income. In the analysis of number of 
prescriptions as a predictor of mortality, we found a 60% 
increase in mortality rate during follow-up for a person 
who, on average, filled one prescription per month 
during the follow-up period compared with a person 
without any prescription fillings (HR=1.60 (95% CI 1.25 
to 2.05)). 



DISCUSSION 
Summary of findings 

In this large Danish twin population we demonstrated 
an inverse social gradient in fillings of prescription 
medicine for all-purpose and system-specific drugs. That 



Madsen M, Andersen PK, GersterM, etal. BMJ Open 201 3;3:e003292. doi:1 0.1 136/bmjopen-201 3-003292 



5 



Open Access 



Table 3 Mean number of prescription fillings after 1 and 
10 years of follow-up in a 5% population-based sample 
and a population of Danish twins, stratified on age and sex 

Type of drugs Nelson-Aalen cumulative hazard 



Age 



5% Population 
(n=137 300) 



Twins (n=24 370) 



Sex 


1 year 


1 0 years 


1 year 


10 years 


All purpose 










30-44 years 










Men 


2.5 


28.7 


2.2 


25.9 


Women 


4.1 


44.5 


3.8 


41.1 


45-59 years 










Men 


3.6 


47.4 


3.3 


43.1 


Women 


6.0 


71.9 


5.7 


68.4 


60+ years 










Men 


6.4 


87.5 


5.8 


82.2 


Women 


8.3 


105.4 


7.2 


95.7 


Circulatory system (ATC-C) 








30-44 years 










Men 


0.1 


3.3 


0.1 


2.9 


Women 


0.2 


4.2 


0.2 


3.6 


45-59 years 










Men 


0.7 


12.6 


0.6 


10.6 


Women 


0.8 


13.8 


0.7 


13.0 


60+ years 










Men 


1.9 


29.6 


1.7 


27.4 


Women 


1.9 


30.0 


1.6 


26.7 


Nervous system 


(ATC-N) 








30-44 years 










Men 


0.8 


9.8 


0.7 


8.6 


Women 


1.1 


14.5 


1.0 


13.3 


45-59 years 










Men 


1.2 


14.6 


1.2 


13.9 


Women 


2.1 


25.0 


2.0 


24.1 


60+ years 










Men 


1.8 


24.4 


1.6 


22.2 


Women 


3.1 


39.2 


2.7 


35.8 


Respiratory system (ATC-R) 








30-44 years 










Men 


0.3 


3.9 


0.3 


3.6 


Women 


0.6 


6.5 


0.5 


5.7 


45-59 years 










Men 


0.5 


5.8 


0.4 


5.0 


Women 


0.8 


9.8 


0.7 


9.1 


60+ years 










Men 


1.1 


14.4 


0.9 


12.5 


Women 


1.2 


14.8 


1.0 


12.8 


Alimentary tract and metabolism (ATC-A) 




30-44 years 










Men 


0.3 


3.5 


0.2 


3.1 


Women 


0.4 


5.0 


0.4 


4.1 


45-59 years 










Men 


0.5 


7.6 


0.5 


6.6 


Women 


0.7 


9.6 


0.7 


8.6 


60+ years 










Men 


1.1 


15.4 


1.1 


15.0 


Women 


1.3 


18.3 


1.0 


16.4 



ATC, Anatomical Therapeutic Chemical classification system. 



is, a greater use of prescription medicine among those 
with a low SEP thus suggesting a worse health status of 
this group. Furthermore, the number of prescription fill- 
ings was shown to be associated with a greater risk of 
mortality. In the intrapair analyses, where shared familial 
factors were controlled per design, the association 
between SEP and prescription fillings attenuated to a 
wide extent within MZ twins and statistically significantly 
less in DZSS twins. For education, there was generally no 
association with prescription fillings in the intrapair ana- 
lyses of MZ twins, except for drugs targeting the nervous 
system, while there was still some association of income 
with all-purpose medication and most system-specific 
drugs. However, no effect of income was seen in the 
intrapair analyses when DDD were used as outcome 
measure. In spite of some variations between system- 
specific drugs, the overall pattern was the same as for all- 
purpose drugs, except for ATC-N drugs and to a lesser 
degree ATC-A drugs, which still showed some effect of 
SEP in the intrapair analyses. 

Interpretation 

The fact that the association was attenuated in the intra- 
pair analyses suggests that underlying factors shared by a 
pair of twins account for some of the observed associ- 
ation between SEP and prescription fillings. However, 
there may be a number of alternative or complementary 
explanations for this finding: first, part of it is likely 
caused by measurement error. Measurement error is 
random and therefore not shared by twins. In the intra- 
pair analysis, variance due to shared factors is removed, 
and a relatively larger proportion of exposure variance is 
thus due to measurement error. This will result in 
increasing attenuation of parameter estimates in the 
intrapair analyses. 33 However, since data come from offi- 
cial registers it can be expected that education and 
income have limited measurement error. Second, the 
behaviour of one twin could be affected by the behav- 
iour of the other twin, that is, 'social interaction' either 
by 'imitation' or 'differentiation'. The implications of a 
potential imitation could be smaller within-pair differ- 
ences, which could be misinterpreted as confounding 
from underlying familial factors. This may be of some 
concern in this study since fillings of prescription medi- 
cine is likely to have a behavioural component. Another 
implication of the discordant twin pair approach is that 
only twin pairs who are discordant on exposure contrib- 
ute with statistical information to the estimation of SEP 
effects in the intrapair analysis. This could potentially 
induce bias in the comparison of the unpaired and 
intrapair analysis, since the two analyses are based on 
only partiy overlapping populations. Since MZ twins are 
less likely to be discordant, the selection bias may be 
most severe in MZ twins. 34 We examined this concern by 
restricting the unpaired analysis only to include groups 
with exposure-discordant twins. Results were similar to 
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Table 4 Ratio (MR) of the mean number of prescription fillings (1995-2005) for all-purpose and system-specific drugs 
according to educational status and income in a 5% population-based sample (137 300) and a population of Danish twins 
(n=24 370), showing results from unpaired and intrapair analyses stratified on zygosity 

Unpaired analysis* Intrapair analysist 

5% Population 



ns 



Twins 

n=137 300 DZSSn=15 788 MZ n=8582 DZSS MZ 



Type of drug MR 95% CI MR 95% CI MR 95% CI MR 95% CI MR 95% CI 

All purpose 
Education 



Secondary 


0.90 


0. 


,90 


to 


0 


.90 


0. 


90 


0 


.86 


to 


0. 


95 


0. 


87 


0. 


81 


to 


0, 


,93 


0. 


88 


0, 


,84 


to 


0. 


,93 


0. 


97 


0 


,91 


to 


1.03 


Tertiary 


0.82 


0. 


,82 


to 


0 


.83 


0. 


85 


0 


.80 


to 


0. 


,91 


0. 


83 


0. 


76 


to 


0, 


,90 


0. 


86 


0, 


,80 


to 


0. 


,92 


0. 


95 


0. 


,86 


to 


1.05 


Income 




































































2nd quartile 


0.78 


0. 


,78 


to 


0 


.78 


0. 


81 


0 


,76 


to 


0. 


,85 


0. 


81 


0. 


.75 


to 


0, 


88 


0. 


86 


0 


,82 


to 


0. 


,91 


0. 


91 


0. 


,85 


to 


0.97 


3rd quartile 


0.72 


0. 


,72 


to 


0 


.73 


0. 


72 


0. 


,69 


to 


0, 


,77 


0. 


74 


0. 


,69 


to 


0. 


80 


0. 


79 


0 


.75 


to 


0. 


,83 


0. 


87 


0. 


,81 


to 


0.94 


4th quartile 


0.70 


0. 


,70 


to 


0 


.70 


0. 


71 


0 


,67 


to 


0, 


,75 


0. 


70 


0. 


,65 


to 


0. 


,77 


0. 


76 


0 


.72 


to 


0, 


,81 


0. 


87 


0. 


,81 


to 


0.94 


Circulatory system 


(ATC-C) 
































































Education 




































































Secondary 


0.95 


0. 


,95 


to 


0. 


.96 


0. 


92 


0 


,85 


to 


0, 


,99 


0. 


87 


0. 


,78 


to 


0, 


96 


0. 


85 


0 


.79 


to 


0, 


,92 


0. 


89 


0. 


,81 


to 


0.97 


Tertiary 


0.79 


0. 


,79 


to 


0 


.80 


0. 


78 


0 


.71 


to 


0, 


87 


0. 


70 


0. 


60 


to 


0, 


,82 


0. 


80 


0, 


.71 


to 


0, 


,90 


0. 


93 


0. 


,79 


to 


1.08 



0.91 


0, 


,90 


to 


0 


.91 


0. 


94 


0 


.87 


to 


1 


.03 


0.94 


0. 


84 


to 


1, 


,06 


1, 


.08 


0. 


.99 


to 


1 


.17 


0 


,95 


0 


,87 


to 


1.05 


0.83 


0, 


,83 


to 


0 


.84 


0. 


86 


0. 


.78 


to 


0 


.94 


0.86 


0. 


76 


to 


0, 


,97 


0 


.97 


0. 


.89 


to 


1 


.06 


0. 


,98 


0, 


,89 


to 


1.08 


0.79 


0, 


,79 


to 


0 


.79 


0. 


82 


0 


.74 


to 


0 


.90 


0.85 


0. 


,74 


to 


0, 


,97 


0, 


.97 


0. 


.88 


to 


1 


.08 


0. 


,94 


0 


,84 


to 


1.05 



Income 

2nd quartile 

3rd quartile 

4th quartile 
Nervous system (ATC-N) 
Education 

Secondary 0.80 0.80 to 0.80 0.87 0.79to0.96 0.77 0.67 to 0.88 0.85 0.77 to 0.94 1.00 0.86 to 1.17 
Tertiary 0.77 0.71 to 0.71 0.77 0.67 to 0.88 0.75 0.61.0.92 0.78 0.67 to 0.90 0.91 0.73 to 1.16 

Income 



2nd quartile 


0.67 0.60 


to 


0. 


.61 


0.67 


0 


,55 


to 


0, 


,68 


0.60 


0. 


,52 


to 


0. 


70 


0. 


68 


0 


.61 


to 


0 


.75 


0. 


75 


0. 


.67 


to 


0.87 


3rd quartile 


0.49 0.49 


to 


0 


.50 


0.47 


0 


.43 


to 


0, 


53 


0.47 


0. 


39 


to 


0, 


,56 


0. 


54 


0. 


.48 


to 


0 


.61 


0. 


59 


0 


,50 


to 


0.70 


4th quartile 


0.44 0.44 


to 


0 


.44 


0.45 


0 


.40 


to 


0, 


51 


0.42 


0. 


,35 


to 


0, 


,50 


0. 


50 


0. 


.45 


to 


0 


.57 


a 


60 


0 


,50 


to 


0.71 


Respiratory system 


(ATC-R) 


























































Education 


























































Secondary 


0.86 0.85 


to 


0 


.86 


0.77 


0 


,67 


to 


0, 


87 


0.82 


0. 


69 


to 


0, 


,98 


0. 


83 


0. 


.72 


to 


0 


.94 


1. 


00 


0. 


,84 


to 


1.21 


Tertiary 


0.S7 0.80 


to 


0 


.81 


0.79 


0 


.68 


to 


0, 


,92 


0.76 


0. 


61 


to 


0, 


,95 


0. 


80 


0. 


.66 


to 


0 


.97 


0. 


77 


0 


,58 


to 


1.03 



Income 



2nd quartile 


0.76 


0.75 to 0.76 


0. 


76 


0 


.66 


to 


0. 


87 


0. 


79 


0. 


,65 


to 


0, 


,98 


0. 


86 


0. 


.75 


to 


0 


.99 


1 


.04 


0 


.86 


to 


1.25 


3rd quartile 


0.68 


0.68 to 0.69 


0. 


74 


0 


,64 


to 


0. 


,86 


0. 


69 


0. 


.57 


to 


0. 


,85 


0. 


85 


0 


.73 


to 


0 


.98 


1 


.02 


0 


.84 


to 


1.24 


4th quartile 


0.63 


0.63 to 0.63 


0. 


69 


0 


,59 


to 


0, 


,81 


0. 


64 


0. 


,52 


to 


0, 


79 


0. 


83 


0 


.71 


to 


0 


.96 


1 


.25 


1 


.04 


to 


1.51 


Alimentary tract and metabolism (ATC-A) 






















































Education 


























































Secondary 


0.84 


0.84 to 0.84 


0. 


85 


0 


,77 


to 


0, 


,94 


0. 


81 


0. 


.71 


to 


0. 


,93 


0. 


84 


0 


.75 


to 


0 


.94 


0 


.98 


0 


.85 


to 


1.14 


Tertiary 


0.68 


0.67 to 0.68 


0. 


72 


0 


.62 


to 


0, 


83 


0. 


64 


0. 


53 


to 


0, 


,77 


0. 


81 


0. 


.68 


to 


0 


.96 


0 


.94 


0 


.76 


to 


1.17 


Income 




























































2nd quartile 


0.74 


0.74 to 0.75 


0. 


73 


0 


.65 


to 


0, 


81 


0. 


83 


0. 


72 


to 


0, 


,96 


0. 


85 


0. 


.76 


to 


0 


.96 


0 


.88 


0 


.77 


to 


1.01 


3rd quartile 


0.64 


0.64 to 0.65 


0. 


63 


0 


.56 


to 


0, 


71 


0. 


77 


0. 


66 


to 


0, 


,90 


0. 


77 


0. 


.68 


to 


0 


.87 


0 


.87 


0 


.75 


to 


1.02 


4th quartile 


0.61 


0.61 to 0.61 


0. 


60 


0 


.53 


to 


0, 


68 


0. 


65 


0. 


,55 


to 


0, 


,77 


0. 


76 


0 


.67 


to 


0 


.88 


0 


.86 


0 


.72 


to 


1.02 



'Unpaired analysis treating twins as individuals while taking interdependence of observations within each individual and within twin pairs into 
account. The interpretation of MR is the ratio of the mean number of prescriptions at a given time of, for example, an individual with a 
secondary or tertiary education compared with a random individual with a primary education. MRs are adjusted for age and sex. Income and 
education are mutually adjusted. 

flntrapair analysis of twins. The interpretation of MR is the ratio of the mean number of prescriptions at a given time of a twin with a 
secondary or tertiary education compared with its co-twin with a primary education. 
Italic typography indicates statistical significance at 5% level. 

ATC, Anatomical Therapeutic Chemical classification system; DZSS, dizygotic same sex twins; MZ, monozygotic twins. 



those of the main analysis in both MZ and DZ twins. 
A further issue relates to the fact that twins who are dis- 
cordant on exposure are special in the sense that they, 
in spite of overlapping genes and rearing environment, 



differ in their SEP. Thus, unshared factors may explain 
these differences. This implies that the intrapair analysis 
may still be confounded by unshared confounding 
factors, and that this confounding may be enforced in 
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Table 5 Results from a linear regression model showing number of daily defined doses (DDD) for all-purpose drugs 
according to educational status and income in a population of Danish twins (n=24 370) stratified on zygosity 

Unpaired analysis* Intrapair analysist 



Educational status 


DZSS n : 


=15 788 




MZ n= 


:8582 




DZSS 






MZ 




Income 


DDD 


95% CI 




DDD 


95% CI 




DDD 


95% CI 




DDD 


95% CI 


Secondary 


-287 


-397 to • 


-177 


-289 


-440 to ■ 


-138 


-245 


-399 to • 


-90 


-143 


-346 to 59 


Tertiary 


-404 


-531 to 


-277 


-424 


-606 to • 


-243 


-344 


-557 to ■ 


-131 


-152 


-440 to 136 


2nd quartile 


-401 


-541 to 


-260 


-275 


-463 to ■ 


-88 


-219 


-383 to ■ 


-54 


18 


-177 to 214 


3rd quartile 


-541 


-680 to 


-402 


-451 


-638 to ■ 


-263 


-367 


-535 to ■ 


-198 


-83 


-287 to 121 


4th quartile 


-505 


-648 to 


-363 


-437 


-630 to ■ 


-243 


-327 


-512 to - 


-142 


16 


-209 to 242 



Results from unpaired and intrapair analyses are shown. 

'Unpaired analysis treating twins as individuals while taking interdependence of observations within each individual and within twin pairs into 
account. The interpretation of, for example, the education estimate is difference in the mean number of DDDs of an individual with a 
secondary or tertiary education compared with a random individual with a primary education. The estimates are adjusted for age and sex. 
Income and education are mutually adjusted. 

flntrapair analysis of twins. The interpretation of for example, the education estimate is the difference in the mean number of DDDs of a twin 
with a secondary or tertiary education compared with its co-twin with a primary education. 
Italic typography indicates statistical significance at 5% level. 
DZSS, dizygotic same sex twins; MZ, monozygotic twins. 



the intrapair analysis.' ' Generally, confounding from 
unshared factors is thought to be more severe when con- 
founders are less correlated than exposure. 33 This may 
or may not be true in our case, since we can only 
hypothesise on the correlation of unmeasured confoun- 
ders. Nevertheless, the correlation of education and 
income (MZ=0.79, DZ=0.58 and MZ=0.49, DZ=0.33) are 
more or less of the same magnitude as phenotypes like 
cognitive ability 35 and body mass index 36 that have previ- 
ously been investigated using similar designs. In any 
case, as shown by Frisell et al S5 in order for increased con- 
founding to explain the attenuation of the association 
observed in this study, the potential non-shared confoun- 
ders would have to create a positive association between 
education and prescription fillings. In most cases, 
however, we would expect a negative association, 
although filling of prescription medicine is indeed a 
complex health outcome and it is theoretically possible 
for a confounder to have a dual effect on prescription 
fillings by (1) increasing the capability of getting the 
relevant prescriptions, for example, due to more regular 
heath checks or better understanding of symptoms (2) 
at the same time reducing the need for prescription fill- 
ings due to better health status. Nevertheless, in a 
Danish context health status does seem to be the main 

• * * 1 0 37 38 

determinant of prescription medicine ' and was 
clearly associated with mortality in this data. Thus, in the 
presence of non-shared confounding we would expect 
the net effect of any confounding to produce a negative 
association between education and prescription fillings 
resulting in a bias away from the null, 33 and therefore it 
could not explain the observed attenuation of effect in 
the intrapair analyses observed in this study. 

The finding of greatest attenuation in MZ twins is 
compatible with the scenario of genetic confounding 
described in the introduction. Since the genetic overlap 
is greater in MZ twins than DZSS twins, one would 
expect a larger attenuation in MZ twins, if genetic 



factors gave rise to the confounding. However, a number 
of other factors could produce a greater attenuation in 
MZ twins: first of all, the aforementioned measurement 
error may be most pronounced in MZ twins, due to the 
fact that increasing correlation in exposure will result in 

r • -33 

a greater attenuation or intrapair parameter estimates. 
In addition, the comparison of strength of association 
between MZ and DZ twins to indicate the source of 
familial confounding relies on the equal environment 
assumption: that is, MZ and DZ twins share their envir- 
onment equally. However, this assumption may be chal- 
lenged, and if it does not hold, that is, if MZ twins not 
only share their genes but also their environment more 
closely than DZ twins, 39 shared environmental factors 
could also produce a greater attenuation of effect in MZ 
twins. However, the interplay between environment and 
genes is complex, and it is conceivable that the environ- 
ment to a certain extent is shaped in response to the 
genetic dispositions of an individual. If, for instance, a 
genetically determined similar appearance of MZ twins 
also elicits a more similar treatment from their peers, 
then these similar environmental stimuli can also be 
ascribed to genetic factors. 40 

Thus, although several factors may contribute to the 
observed attenuation of association between SEP and 
filling of prescription drugs in the intrapair analyses and 
to the differences in strength of association between MZ 
and DZ twins, the findings suggest that a part of the 
observed social inequality in prescription fillings is 
explained by either genetic or environmental factors 
shared by twins. 

The finding of a persistent association of particularly 
income with fillings of prescriptions for ATC-N drugs 
clearly stood out from the general results. We therefore 
carried out a post hoc subanalysis of specific categories 
of ATC-N drugs to test whether it were certain types of 
drugs driving the association. The analysis revealed that 
it was mainly for analgesics and psycholeptics that the 
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association persisted, while the association was somewhat 
attenuated for psychoanaleptics (primarily antidepres- 
sants). This finding may support the notion of reverse 
causation, that is, that health status influences SEP 
rather than the other way around, since it is likely that 
diagnoses such as chronic pain disorders, psychoses, 
anxiety and bipolar disorders associated with these types 
of drugs may affect the work ability and thus the income 
of a person. Income is much more sensitive to reverse 
causation than education, which is usually obtained 
early in life and remains constant once it is obtained. 

Existing literature 

In this study we argue that fillings of all-purpose and 
system-specific drug prescriptions to some extent serve 
as an indicator of general health status and of conditions 
such as ischaemic heart disease, hypertension, COPD 
and diabetes. The finding of weaker effects of SEP, when 
underlying familial factors were taken into account, 
accords with a large part of the life-course literature that 
supports the importance of early life factors in the aeti- 
ology of chronic disease in adulthood. 41 42 

In addition, a number of previous twin and sibling 
studies have investigated the association between various 
SEP indicators and health outcomes. While some studies 
have shown strong or some evidence of a causal effect of 
adult SEP 43-47 other suggest no or little effect of SEP, 
when familial factors have been adjusted for. 22 48-51 
Other studies show mixed results. 21 23-25 47 52-54 
However, the studies vary with respect to a number of 
important parameters hampering overall comparisons. 
Most importantly, investigated health outcomes differ 
and may have different aetiologies and show different 
social gradients and patterns of familial confounding, 
but even for specific-health outcomes there is no clear 
picture of evidence. For example, for mortality, two 
sibling 45 46 and two twin studies 21 44 have investigated 
the association with education. Although all studies 
support some effect of education after controlling for 
shared family factors, important differences exist: the 
twin study by Madsen et al 21 demonstrated an effect of 
education within twins only in older cohorts of men, but 
not in younger men and not in women. The influence 
of shared family factors was thus suggested in the latter, 
but not in the first subgroup of this study. The two 
sibling studies also suggested some effect of shared 
family factors, while the twin study by Lundborg 43 did 
not support any such effect. For self-rated health, the 
majority of the existing twin studies 24 25 3 have shown 
mixed results, although Osier et at 7 and Lundborg 43 did 
find a statistically significant association between social 
class respective education and self-rated health in the 
intrapair analyses. One important reason for the incon- 
sistencies may be that many studies have suffered from 
power limitations in the intrapair analyses due to a small 
number of events. 21 23 25 47 53 This is not a concern in 
this study, where almost all of the exposure discordant 
twins contributed with information to the intrapair 



analyses, thus providing some of the most powerful twin 
results on social inequality in health to date. 

Strengths and limitations 

The twin approach is clearly a strength of this study, 
since it makes it possible to adjust for unmeasured back- 
ground factors shared by twins that may underlie the 
association between SEP and health. Furthermore, the 
large sample size and the long follow-up allow a great 
number of events to be analysed. In addition, the use of 
population-based registers has the advantage of relatively 
little selection and measurement error in the study 
population. Finally, this study not only focuses on pre- 
scription medicine in general, but also on drugs target- 
ing specific organ systems which may give us a better 
idea of potential mechanisms. 

Among study limitations is the fact that filling of pre- 
scription medicine is not a perfect measure of health 
status. We do, however, argue that the observed associ- 
ation with mortality and the inverse relationship with 
education and income likely reflect a greater need for 
medicine among people of lower SEP. In case that access 
to healthcare also plays a role, presumably favouring 
people of high SEP, the social gradient in health status 
in this study could be underestimated. Although data 
did not contain information on over-the-counter medica- 
tion, it only constitutes a minor part of the total drug 
consumption in Denmark, 7 and since we were interested 
in drugs as an indicator of health status, it could be 
argued that the most valid approach is to focus on drugs 
that require a formal medical indication (ie, 
prescription). 

As models for counting processes only allow one event 
at a time, the number of prescription fillings was 
defined as the main outcome of interest. Therefore, it 
was reassuring that the subanalysis on DDD led to 
similar results. 

Although the distribution of SEP indicators and pre- 
scription fillings were not entirely identical in twins and 
the population-based sample, the associations between 
SEP and prescription fillings were indeed, which sup- 
ports generalisability to the general population in line 
with previous studies that show that health and mortality 
in twins are similar to the background population. 55 56 
Furthermore, the fact that the analysis based only on dis- 
cordant twin pairs show marginal results similar to those 
based on the total twin sample was reassuring. 

Conclusion and perspectives 

In this study, based on a comprehensive twin cohort, we 
observed a social inequality in fillings of prescription 
medicine, showing fewer prescription fillings among 
people with high SEP. This association attenuated within 
twin pairs for most types of prescription fillings, except 
those targeting the nervous system. The attenuation was 
most pronounced in MZ twins. Measurement error can 
contribute to the observed attenuation of association, 
but it is likely to account only for a minor fraction, 



Madsen M, Andersen PK, GersterM, etal. BMJ Open 201 3;3:e003292. doi:1 0.1 1 36/bmjopen-201 3-003292 



9 



Open Access 



considering the exposure and outcome measures in this 
study. As for potential non-shared confounding we 
would, if anything, expect it to bias the estimates away 
from the null. 

We therefore interpret the findings as an indication 
that at least a part of the observed social inequality in 
prescription fillings is explained by shared familial 
factors. The implications of these findings are by no 
means a disqualification of SEP as a determinant of 
health, but rather an indication that family matters in 
the production of social health inequalities; most likely 
through early environment, and possibly also genetic 
factors. Recognising biological components and early 
life exposures as important explanatory factors may 
improve our understanding of the multifaceted aetiology 
of health inequalities over the life course and increase 
the potential for effective and targeted interventions. 
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